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The accurate measurement of extragalactic distances is a central challenge 
of modern astronomy, being required for any realistic description of the age, 
Q\ ', geometry and fate of the Universe. The measurement of relative extragalactic 
distances has become fairly routine, but estimates of absolute distances are 
rare0. In the vicinity of the Sun, direct geometric techniques for obtaining 
absolute distances, such as orbital parallax, are feasible, but heretofore 
such techniques have been difficult to apply to other galaxies. As a result, 
uncertainties in the expansion rate and age of the Universe are dominated b 
uncertainties in the absolute calibration of the extragalactic distance ladde 
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Here we report a geometric distance to the galaxy NGC4258, which we infer 



an important role in future distance-scale calibrations. 



from the direct measurement of orbital motions in a disk of gas surrounding 
the nucleus of this galaxy. The distance so determined - 7.2 ± 0.3 Mpc - is the 
most precise absolute extragalactic distance yet measured, and is likely to play 

in 

NGC4258 is one of 22 nearby AGN known to possess nuclear water masers (the 
microwave equivalent of lasers). The enormous surface brightnesses 10 12 K), small 
sizes (<z 10 14 cm), and narrow linewidths (a few km s -1 ) of these masers make them ideal 
probes of the structure and dynamics of the molecular gas in which they reside. VLBI 
observations of the NGC4258 maser have provided the first direct images of an AGN 
accretion disk, revealing a thin, subparsec-scale, differentially rotating warped disk in the 
nucleus of this relatively weak Seyfert 2 AGN Iffi!. Two distinct populations of masers 
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exist in NGC4258. The high-velocity masers amplify their own spontaneous emission and 
are offset ±1000 km s -1 and 4.7-5.1 mas (0.16-0.28 pc for a distance of 7.2 Mpc) on either 
side of the disk center. The beautiful Keplerian rotation curve traced by these masers 
requires a central binding mass (M), presumably in the form of a supermassive black hole, 
of (3.9 ± 0.1) x 10 7 (£>/7.2 Mpc)(sinz s /sin82)- 2 solar masses (M ) where D is the distance 
to NGC4258 and i s is the disk inclination. Because the high-velocity masers lie in the plane 
of the sky, they should to first order remain stationary as the disk rotates. The systemic 
masers, on the other hand, are positioned along the near edge of the disk and amplify 
the background jet emission evident in Figure 1 111. A fundamental and as-yet untested 
prediction of the maser disk model is that the systemic masers should drift with respect to 
a fixed point on the sky by a few 10 /zasyr -1 as the disk rotates at ~ 1000 km s^ 1 . 




Fig. 1. — See captions at end of text. 

NGC4258 was observed at 4-8 month intervals between 1994 and 1997 with the VLBA 
of the NRAO in order to search for the expected motions. Since the maser emission is 
essentially continuous across the envelope of systemic maser emission, we are forced to 




Fig. 2. — See captions at end of text. 
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rely on structure in the systemic spectrum to isolate and track individual maser features. 
The assumption that distinct peaks in the spectrum correspond to individual clumps of 
gas is justified by the successful tracking of maser accelerations in single-dish monitoring 
programs, and leads to an identification of 20-35 potentially trackable systemic masers in 
each epoch. These features are too densely packed in position and velocity to be individually 
tracked with any reliability. However, the resolution and sampling are sufficient to detect 
bulk rotation in the system, and we have developed a Bayesian pattern-matching analysis 
tool to track inter-epoch shifts in the positions and velocities of the systemic masers, as a 
whole!. The analysis assumes the systemic masers are randomly and narrowly scattered 
about an average radius, (r s ), of 3.9 mas, as indicated by the global disk-fitting analysis 
(see Figure 1). The precise magnitude of the radial scatter is set so as to maximize the 
overall likelihood of the tracking analysis. In order to evaluate the likelihood of a given 
bulk proper motion ((9 X )) or acceleration ((vi os )), the pattern-matching procedure must 
compute the likelihood that each individual maser has in fact moved by (6 X ) or (vi os ). 
This leads to robust estimates for the "trackability" of each maser. Figure 2 shows the 
best-fitting acceleration and proper motion tracks for the most reliably trackable systemic 
masers. Figure 3 shows the overall probability density functions (PDFs) for (vi os ) and (0 X ). 
The PDFs indicate a highly significant detection of bulk motion in the disk and from them 
we conclude (vi os ) = 9.3 ± 0.3kms _1 yr -1 and (8 X ) = 31.5 ± 1 //asyr -1 , where these and all 
subsequent uncertainties are la values. The latter result is consistent with expectations 
and is the first detection of transverse motion in the NGC4258 accretion disk. We note that 
the pattern-matching algorithm has been verified on a number of simulated datasets with 
feature densities and spectral and spatial resolutions comparable to those of the true data. 

In order to convert the maser proper motions and accelerations into a geometric 
distance, we express (6 X ) and (vi os ) in terms of the distance and four disk parameters: 
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Here _D 6 is the distance in Mpc, a s is the disk position angle (East of North) at (r s ), and 

is M/D sin 2 i s as derived from the high- velocity rotation curve and evaluated at 
D = 7.2 Mpc and i s = 82.3° (in units of 1O 7 M ). Q s = (GM 7 . 2 / (r s ) 3 ) 1/2 is the projected 
disk angular velocity at (r s ) as determined by the slope of the systemic position-velocity 
gradient (in units of kms -1 mas -1 ; see Figure 1). a priori estimates for each of these disk 
parameters, derived directly from the positions and velocities of the masers, are included in 
the denominators of each of the terms of equations p] and |2|. 
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Fig. 3. — See captions at end of text. 
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When the a priori disk parameter estimates are used, the proper motions and 
accelerations yield independent distance estimates, through equations [i] and |2], of 
7.2 ± 0.2 Mpc and 7.1 ± 0.2 Mpc, respectively. The quoted uncertainties are effectively the 
uncertainties in (9 X ) and (vi os ) recast in terms of distance, and as such are purely statistical 
in nature. The excellent agreement between the proper motion and acceleration distances for 
a priori values of the disk parameters is an impressive confirmation of the a priori disk model 
itself, and establishes the NGC4258 Keplerian disk as a fully self-consistent, dynamical 
model incorporating the positions, LOS velocities, proper motions, and accelerations of all 
the NGC4258 masers. We note that a preliminary geometric distance estimate, based on 
accelerations alone and a single VLBA epoch, yielded D = 6.4 ± 0.9 Mpc §. Hence, the old 
and new distances are consistent at the la level. The discrepancy between the estimates 
is ultimately explained by the fact that the original distance estimate assumed an average 
systemic maser radius about 10% larger than the value indicated by the newer data and the 
more sophisticated disk models. 

Uncertainties in the disk parameters contribute to systematic uncertainties in the 
distance estimate. We derive a composite geometric distance estimate using equations [I] 
and |2|, and from the (vi os ) and (6 X ) PDFs of Figure 3 and the a priori estimates for the 
disk parameters and their associated uncertainties. The result is a geometric distance 
estimate of 7.2 ± 0.3 Mpc, where the quoted uncertainty now incorporates all statistical 
terms associated with tracking motions in the disk as well as systematics arising from disk 
parameter uncertainties. A ~ 5% uncertainty in Q s is the dominant contributor to the 
latter, producing fractional uncertainties in the acceleration and proper motion distances 
of 6.7% and 1.7%, respectively. In total, disk- model systematics contribute an additional 
0.26 Mpc (in quadrature) to the distance error budget derived from purely statistical 
considerations. 

The NGC4528 geometric distance is the most precise, absolute extragalactic distance 
measured to date and, being independent of all other distance indicators, it represents 
an important new calibration point for the extragalactic distance ladder. The geometric 
distance is consistent with pre-existing H-band Tully-Fisher (7.1 ± 1.1 Mpc blue 
Tully-Fisher (7.9 ± 1.8 Mpc 0), and luminosity class (8.4 ± 2.2 Mpc §) distance estimates, 
all of which rely on the Cepheid Period-Luminosity relationship for absolute calibration. 
Most importantly, efforts are underway to determine directly a Cepheid distance to 
NGC4258 from Hubble Space Telescope (HST) observations of the galaxy. The NGC4528 
geometric distance is presently the most precise means for directly calibrating distances 
obtained by HST Cepheid observations. 

We emphasize that the above error budget considers statistical and systematic 
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uncertainties within the framework of a thin Keplerian disk in which the masers trace the 
orbital motions of discrete clumps of gas. As always, it is difficult to estimate any additional 
systematic uncertainties that might exist as a result of imperfections in the model itself. 
The potential impact of eccentricity on the distance error budget depends on the assumed 
distribution of accretion disk eccentricities in AGN in general. Viscous dissipation within 
such disks is expected to circularize orbits on relatively short timescales, and detailed 
modeling of the optical emission lines from a large sample of AGN suggests eccentricities 
less than about 0.5 111. The eccentricity of the NGC4258 disk is further constrained by 
the symmetry of the maser emission about the disk center in both position and velocitylll. 
These additional constraints lead to an expected eccentricity of zero with a probable error of 
0.1, a negligible bias in the distance estimate, and a systematic uncertainty in the distance 
of 0.4 Mpc. Hence our distance estimate, including this uncertainty in the eccentricity, 
is 7.2 ± 0.5 Mpc. Finally, we cannot unambiguously rule out contamination of the maser 
dynamics by some non-kinematical contribution, such as traveling density waves within the 
disk. However, given the complexity (30 systemic masers across 8° of disk azimuth) and the 
stability 70% of the features persisting) of the pattern we have tracked, orbital motion 
is certainly the simplest explanation. 
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CAPTIONS 

Figure 1. - The upper panel shows the best-fitting warped disk model superposed on 
actual maser positions as measured by the VLBA of the NRAO, with top as North. The 
filled square marks the center of the disk, as determined from a global disk-fitting analysis!. 
The filled triangles show the positions of the high-velocity masers, so called because they 
occur at frequencies corresponding to Doppler shifts of ~ ±1000 km s" 1 with respect to 
the galaxy systemic velocity of ~ 470 km s" 1 . This is apparent in the VLBA total power 
spectrum displayed in the lower panel. The inset shows line-of-sight (LOS) velocity versus 
impact parameter for the best-fitting Keplerian disk, with the maser data superposed. 
The high- velocity masers trace a Keplerian curve to better than 1%. Monitoring of these 
features indicates that they drift by less than ~ lkms^yr" 1 iiH and requires that they 
lie within 5-10° of the midline, the intersection of the disk with the plane of the sky. The 
LOS velocities of the systemic masers are centered about the systemic velocity of the galaxy. 
The positions (filled circles of upper panel) and LOS velocities of these masers imply they 
subtend about 8° of disk azimuth centered about the LOS to the central mass, and the 
observed 8-10 km s -1 yr" 1 acceleration of these features!! unambiguously places them along 
the near edge of the disk. The approximately linear relationship between systemic maser 
impact parameter and LOS velocity demonstrates that the disk is very thinl^l (aspect ratio 
<z 0.2%) and that these masers are confined to a narrow annulus in the disk. The magnitude 
of the velocity gradient (Q s ) implies a mean systemic radius, (r s ), of 3.9 mas which, together 
with the positions of the high- velocity masers, constrains the disk inclination, i s , to be 
~ 82 ± 1° (90° for edge-on). Finally, VLBA continuum imagesBli! are included as contours 
in the upper panel. The 22-GHz radio emission traces a sub-parsec-scale jet elongated along 
the rotation axis of the disk and well-aligned with a luminous, kpc-scale jetS. 

Figure 2. - Line-of-sight (LOS) velocities (a) and right ascensions (b) at the peaks 
of the systemic maser spectrum for each of the five VLBA epochs. Only those features 
deemed reliably trackable by the pattern-matching analysis are included. All epochs 
included the Very Large Array, phased to act as a single large aperture. In addition, the 
final epoch utilized the Effelsberg 100-meter telescope. In each epoch the VLBA correlator 
provided cross-power spectra with channel spacings of 0.22 km s" 1 . The NGC4258 masers 
are characterized by linewidths of ~ 1 — 3kms _1 . All spectra were phase-referenced, 
via self-calibration, to a strong systemic maser to stabilize the interferometer against 
atmospheric pathlength fluctuations, and synthesis images were constructed for each 
spectral channel using conventional restoration techniques^. Error bars have been foregone 
in order to avoid clutter. The uncertainties in individual LOS velocity estimates are 
dominated by line blending in the spectrum. There is an average scatter of about 0.4kms _1 
about the best-fitting acceleration tracks. The masers are spatially unresolved and relative 
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positions have been measured to a precision of ~ 0.56 B /SNR, where SNR is the signal to 
noise ratio and 0# represents the 0.6 x 0.9 mas, approximately North-South synthesized 
beam. Relative positional accuracies typically ranged from 0.5 to 10 //as. All positions 
are relative to a fixed point along the systemic position-velocity gradient. Our ability to 
precisely align this structure amongst all epochs suggests that it does indeed remain fixed in 
time. The best-fitting acceleration and proper motion tracks (solid lines) indicate average 
drifts of 9.3kms _1 yr _1 and 31.5;uasyr _1 in velocity and position, respectively. The scatter 
in the individual proper motions and accelerations about these average values is consistent 
with a 0.2 mas scatter in the radii of the systemic masers about (r s ). 

Figure 3. - Systemic maser bulk proper motion ((0 X ); a) and acceleration ((vi os )] 
b) probability density functions (PDFs) as derived using the Bayesian pattern-matching 
analysis described in the text. The curves were generated using all the maser features 
in each epoch. The uncertainties in (9 X ) and (vi os ) as derived from these PDFs include 
measurement uncertainties in the maser positions and velocities, but they are dominated by 
ambiguities in tracking specific maser features. 



